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This work addresses the question of whether low-lying excitations in conjugated polymers are com- 
prised of free charge-carriers or excitons. States are characterised as bound or unbound according 
to the scaling of the average particle-hole separation with system size. We critically examine other 
criteria commonly used to characterise states. The polymer is described by an extended Hubbard 
model with alternating transfer integrals. The model is solved by exact diagonalisation and the 
density matrix renormalisation group (DMRG) method. We demonstrate that the DMRG accu- 
rately determines excitation energies, transition dipole moments and particle-hole separations of a 
number of dipole forbidden {Ag) and dipole allowed (-Bu) states. 

Within a parameter regime considered reasonable for polymers such as polyacetylene, it is found 
that the charge gap, often used to define the exciton binding energy, is not a good criterion by 
which to decide whether a state is bound or unbound. The essential non-linear optical state mAg 
is found to mark the onset of unbound excitations in the Ag symmetry sector. In the Bu symmetry 
sector, on the other hand, it is found that all low lying states are unbound and that there is no well 
defined uBu state. That is, the IB^ state marks the onset of unbound excitations in this sector. 



PACS numbers: 71.10.F, 71.20.R, 71.35 



I. INTRODUCTION 

The current interest in conjugated polymers lies to 
a large extent in their optical propertiesEl. Conjugated 
polymers exhibit strong luminescence, and large and ul- 
trafast nonlinear optical (NLO) responseH. This has lead 
to technological opportunities and, for instance_polymer 
light-emitting diodes are now widely produceda. These 
optical properties are associated with the delocalised tt- 
electron system of the conjugated polymers and, in par- 
ticular, the low-lying excitations. However, the nature of 
these excitations is not fully understood, and has been a 
subject for fundamental research on conjugated polymers 
in recent yearsEl. 

A central issue is whether low-lying excitations are 
comprised of free charge-carriers or excitons. If the 
Coulomb interaction between the oppositely charged par- 
ticles and holes is strong, excitons are formed, i.e. bound 
particle-hole pairs, in which the motions of the particle 
and the hole are strongly correlated. On the other hand, 
if the Coulomb interaction is effectively screened, then 
the particles and holes are only very weakly bound and 
move essentially independently as free charge-carriers. 

Different criteria have been used in the literature to 
discriminate between the states, leading to contradic- 
tory conclusions concerning the nature of the low-lying 
states- Some commonly used criteria are: the charge \ 
iJgQO'El; and the essential NLO states, mAg and n_B„ 



The charge gap is defined as the sum of the energies for 
removing and adding an electron to the neutral system: 



Eg = E{N + 1) + E{N - 1) - 2E(N). 



(1) 



where E{N) refers to the ground state energy of the 
neutral system (where the band is half filled). For the 
Hubbard model with only on-site Coulomb interaction, 
the charge gap coincides with the lowest optical exci- 
tation, i.e. the optical gapQ. This is, of course, also 
true for the tight-binding (or Hiickel) model, which is an 
independent-electron model, where the charge gap is the 
onset of the delocalised states in the conduction band. 
However, this is generally not true in the case of longer 
range Coulomb interactions since there are states below 
the charge gap. It has been argued that these states 
are exciton states since they, for instance, appear ener- 
getically within the tight-binding band gap. Thus, the 
charge gap has been used to discriminate between the 
states: states below the charge gap correspond to exci- 
tons and states above the charge gap correspond to free 
charge-carriers. Consequently, for a state of energy E, 
the binding energy E\^ is defined as 



Eh = E„ 



E. 



(2) 



We believe that this criterion can be seriously criticised: 
(i) the charge-gap energy is in general not an eigenenergy 
of the system, but rather mixes ground state energies of 
three, all differently charged, systems; (ii) the criterion 
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does not distinguish between different symmetry sectors; 
and (iii) the criterion does not directly measure the mo- 
tion of particles and holes, but is based on total energies 
only. In this work, we have directly calculated the rela- 
tive motion of particles and holes for different states in 
each symmetry sector in order to identify excitons. 

It should be noted that the use of the Hartree-Fock 
bandgap in the literatureEl is merely an approximate way 
of calculating the charge gap. This follows directly from 
Koopmans' theorem, and we will therefore not treat the 
Hartree-Fock bandgap as separate criterion for excitons. 
It is well known that the Hartree-Fock bandgap system- 
atically overestimates the true charge gap. 

In works on NLO properties of conjugated polymers 
the most kpapttapt channels for such processes have been 
identified@ll3'tiHl3, leading to a phenomenological model 
for third-order nonliiiearity that is based on only the four 
most essential statesQ. These states are the ground state 
lAg, the lowest dipole allowed state IB^, the mAg and 
the nBu- The roAg is defined as the state which has the 
strongest dipolellj coupling (or transition moment) to the 
IBu, and the nB^ is defined as the state which has the 
strongest dipole coupling to the mAg, apart from the 

In addition, it was found that there is sudden increase 
in the particle-hole separation at the mAg and the nBu, 
i.e. all states belowJhese states have more tightly bound 
particles and holesEl It has therefore been argued that 
the mAg and the nB^ are the lowest lying free charge- 
carrier states and useful criteria for the identification of 
excitons. The binding energy is defined analogously to 
(^, but with the charge gap replaced with the mAg or 
nBu energy depending on which symmetry sector is of 
interest. 

Although intriguing, it is not clear that states below 
the mAg and the nBu really correspond to excitons. The 
smaller particle-hole separation may simply be a conse- 
quence of system confinement. In principle one needs to 
study an infiaite system to resolve this issue, whereas 
the results oiul pertain to an oligomer of iV = 8 carbon 
atoms. Another approach is to look at how the particle- 
hole separation scales with the system size N. In this 
work, we have calculated the particle-hole separation for 
different system sizes. The particle-hole separation scales 
differently for free charge-carriers and bound states. 

In another study, the particle- hale separation was stud- 
ied as a function of system siza3, with electron corre- 
lation treated within the singly excited configuration- 
interaction (SCI) approximation. It was found that 
states of energy significantly below (above) the Hartree- 
Fock band gap were bound (unbound). Excitons are 
many body excitations and, for instance, it has been 
found that the exciton binding energy is sensitive to the 
inclusion ofJiigher- order correlation through perturba- 
tion theorjij'lla. It is therefore important to treat the 
electron correlation accurately when assessing exciton 
criteria. In this work, we have used exact diagonalisation 
of the Hamiltonian for systems {N < 10), and the Den- 



sity Matrix Renormalisation Group method for longer 
systems {N < 50). 

The methodology of this work is outlined in Sec. ||, 
which contains definitions of calculated quantities and a 
description of the computatio nal m ethods. The model 
Hamiltonian is defined in Sec. II A , particle- hole separa- 
tion definitions are given in Sec. [IB where the basis of 
the scaling analysis is outlined, the ionicity is defined in 
Sec. H C and t he nu merical solution of the model is de- 
scribed in Sec. HD. A more detailed derivation of the 



scaling behaviour in the non-interacting limit is given in 
Appendix A. Results are presented and discussed in Sec. 
III. Th e acc uracy of the DMRG method is demonstrated 
in Sec. HI A, the energy spectra and transition moments 
are given in Sec. IIIB , results for the particle-hole sepa- 
ration follow in Sec. Ill C , whilst ionicity results are given 



in Sec. Ill D . Finally, conclusions are given in Sec. 



[V 



II. METHODOLOGY 



A. Model 



As a generic model for conjugated polymers we study 
the extended Hubbard model with alternating hopping 
integrals and on-site and nearest-neighbour electron- 
electron interactions: 

H ^ -tJ2[^ - (-l)^<5](cLc.+i. + h.c.) + 

I - l){n, - 1) + VY,ifl^ ~ 1 - 1), (3) 

i i 

where Cia annihilates a 7r-electron with spin a on site i 
and "hi is the occupation number operator for site i. We 
have studied systems with an even number of sites N, 
and open boundary conditions. The Hamiltonian is a 
paradigm for conjugated polymers in which excitons can 
exist. 

For convenience we take t = I, which sets the energy 
scale. Values for the dimerisation in the range 0.07-0.15 
have been proposed for polyacetYlene, polydiacetylene, 
and poly(para-phenylenevinylene)ll3'cl. We choose S — 
0.1 as a typical value for conjugated polymeps. Optical 
absorption data suggest a U value of 2.25-2.75liIl, whereas 
ab Jnitio calculations suggest a slightly higher value of 
3.6113. We have taken U = 3 as a reasonably realistic 
value. For the nearest nedshbour interaction we use the 
standard value V — OAUuu. Thus, the chosen values of 
the parameters are: t = 1, 6 = 0.1, U = 3, and V = 
1.2. These values are used throughout the paper, unless 
otherwise stated. 

Besides conserving the total number of particles, 
the Hamiltonian possesses spatial symmetry, C2h G 
{Ag,Ba}, charge conjugation symmetry. J G {+,—}, 
and spin symmetry, e.g. S G {0,1,2.. .}Ej. In this pa- 
per we are concerned with neutral, singlet states (with a 
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half-filled band and S = 0), except for the charged states 
used in determining Eg. The singlet ground state is even 
under spatial inversion and charge conjugation, and is 
denoted l^A^, or simply lAg. For this work we only 
need to consider the ground state symmetry sector and 
the sector to which it is dipole coupled. Since the dipole 
operator is odd under inversion and charge conjugation, 
states in the dipole coupled sector are denoted j^B^ , or 
simply j-Bjj, where j is the state number. 



B. Particle- hole separation 

In this paper, we consider the relative motion of par- 
ticles and holes as a direct way of identifying excitons. 
In particular, we calculate the average particle-hole sep- 
aration. Excitons have small particle-hole separations 
which remain finite as the system size is increased. By 
contrast, the average separation between two free charge- 
carriers increases indefinitely with system size. Indeed, 
for a completely independent particle and hole, each with 
a uniform probability (l/N) of being on any given site, 
the leading term in the average separation is proportional 
to N. p 

In referenced Guo et al. used the density-density cor- 
relation function C{i,j) as a signifier which distinguishes 
bound and unbound states: 



C{hj) = {{fit - - {fi]))) 



(4) 



C{i,j) correlates a charge fluctuation on site i to a charge 
fluctuation on site j. A positive value means that an ex- 
cess (deflcit) of charge on site i correlates with an excess 
(deficit) on site j. A negative value correlates an excess 
with a deficit. In the context of excitons, the concept of 
quasi-particles is a way of representing an aggregation of 
electrons which leads to an excess of charge in a region, 
which may extend over several sites. Likewise, a hole 
represents deficit of charge. Now, for sufficiently large 
N, C{i,j) is negative (positive) for \i ~ j\ odd (even). 
We thus consider the odd distance contributions when 
defining the probability distribution from which to mea- 
sure the separation of particles and holes. Other defini- 
tions, utilising both positive and negative contributions, 
or |C(i,j)| are possible, and yield the same qualitative 
results. 

InQ it was found that there is a noticeable change in the 
nature of the decay of C{i,j) at the mAg and the 
particle-hole separation being considerably larger than 
for the lower lying states. However, calculations were 
restricted to = 8. We have used the DMRG method 
to study substantially larger systems (up to A^ = 50). 
Most importantly, this makes it possible to discriminate 
between bound and unbound states from the scaling of 
C{i,j) or the particle-hcde separation with A^. 

We define an averagedtj, centered correlation function 
for a system of size N with open boundary conditions as: 



Cn{1) 



r<(N~i N+i \ 

1 \fi( N-l+l N+l+l ' 

2 L^V 2 ' 2 



I even 

C(iv^,i^)] I odd. 

(5) 



As discussed, we define the average particle-hole separa- 
tion (in units of chemical bonds) by regarding the nega- 
tive values of Cat (^) as a probability distribution viz. 



_j:l~o'ii\CN{i)\-c^{i)) 



Ef=o\\CN{i)\ 



Cjv(O) 



(6) 
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FIG. 1. Particle-hole separation 
IBu states in the non-interacting limit {U — V 
as a function of 1/A*'. 



of the lAg, 2Ag and 



0) plotted 



It is instructive to consider the two extreme cases of 
very weak and very strong electron-electron interaction, 
since the identification of bound and unbound states is 
unambiguous in these cases. In the non-interacting limit 
{U = V = 0), we combine exact analytical results with 
exact diagonalisations of systems of up to 4098 sites in 
order to determine the scaling of the correlation function 
and particle-hole separation with A^. Details are given in 
Appendix A. The main point is that, as can be seen from 
Fig. |l|, the ground state {lAg) and the unbound states 
{IBu and 2Ag) have particle- hole separations which scale 
very differently with A^. For the ground state we have 
exponential convergence 



(O^ = «[l + 0(e-"^)] 



(7) 



and for the unbound excited states we have very slow 
convergence 



{l)^^b[l + 0{N~^)] 



(8) 



In the limit of stroag interactions, {U, V ^ t), it was 
shown by Guo et alB that the states can be identified 
solely from their total energy. At t = 0, the anal- 
ysis is trivial: A particle (hole) is a doubly-occupied 
(empty) site. The ground state has zero energy, bound 
states occur a.t U — V, and unbound states at U. Since 
the Hamiltonian consists only of the U-V potential, the 
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bound states have a particle and a hole next to one an- 
other, the correlation function vanishes for I > 1, and 
the particle-hole separation converges immediately with 
N, i.e. = 1. For unbound states, at energy [/, it is 
easy to show that increases linearly with N, i.e. as 
for independent particles and holes. 

In addition to the correlation function, we consider a 
reduced correlation function: the difference between the 
correlation function of an excitation, Cn(1), and that of 
the ground state, c'j^^\l): 



(GS), 



(9) 



The motivation for studying C^'''*^ (/) is that it only mea- 
sures changes in the charge fluctuations induced by an 
excitation. That is, the effect of creating particles and 
holes, whose motion we are interested in. By analogy 
with (||), we also define a reduced particle-hole separa- 
tion 



C. lonicity 

In the limit of strong electron-electron interaction, 
U ^ t, the ground state is a pure spin-density- wave, 
i.e. it is solely a linear combination of covalent configu- 
rations (or Slater determinants), where there is exactly 
one electron on each site. It is therefore natural, in this 
limit, to define a particle as a doubly occupied site, and 
a hole as an unoccupied site. Configurations which have 
particles and holes are referred to as ionic. A general 
state can ba-written as * = ^f(™vaicnt) _^ ^ (ionic) ^ 

Following^, we_define the ionicity as the number of 
particle-hole pairsEiJ. Thus, single particle-hole pair ex- 
citons are constructed from singly ionic configurations, 
biexcitons are constructed from doubly ionic configura- 
tions etc. The ionicity of a state is defined as the expec- 
tation value of the ionicity operator /: 



I) 



(10) 




FIG. 2. Reduced particle-hole separation (l)^'^'''^'' of the IBu 
and 2Ag states in the non- interacting limit {U = V — 0) 
plotted as a function of A'^. 

In the non-interacting case, the reduced particle-hole 
separation for unbound states scales linearly with N, 
as shown in Fig. |^. As mentioned above, this is how 
the separation of a completely independent particle and 
hole scale. In the strongly correlated limit, C]^ (/) = 
and so {l)'"^'^'^^ = (Oat and hence (Oiv'^^'' = 1 for bound 

states and (Oiv^^'' ^ ^'^^ unbound states. That is, un- 
bound states scale in the same way in both limits. This 
suggests that, for general interaction strength, the re- 
duced particle-hole separation scales linearly for unbound 
states. 

In summary, for unbound states, the particle-hole sep- 
aration converges slowly with N and the reduced 

particle-hole separation {1)%'^'^^ diverges linearly. By con- 
trast, for bound states, both (Z)^ and (Ojv'^'^'' converge 
rapidly to a finite value. 



the charge gap-and IBu 
polarisabilitie£3, and os- 



In addition to the ionicity, we have calculated the aver- 
age ionic particle-hole separation (Z)^°'"'^'' by considering 
only the ionic part of the wavefunction. In contrast to 

(Oat: (OaT""^'' allows us to probe changes within \I>('°""=) 
alone, rather than at the same time taking into account a 
change in the relative weig htS of ^-(covalent) ^^^^ ^(ionic)^ 



D. Computational methods 

Equation (||) can be comfortably solved by exact di- 
agonalisation for systems of up to iV = 12 sites. For 
longer chains we turn to the Density Matrix Renormali- 
sation Group (DMRG) methodB. The DMRG has been 
applied to (^) in calculatioiis^of ground state and triplet 
energies charge dcnsitiesu; 
energiesQcJ, the 2-4^ energyE 
cillator strengthsE^. 

In thiS|_work we apply an infinite lattice DMRG 
algorithmEj to find the charge gap and a number of Ag 
and Bu states of (^). In addition to excitation energies, 
we find the transition dipole moments between the Ag 
and Bu states as well as the density-density correlation 
function (^) and hence the particle-hole separation (||) 

and reduced separation (Oiv'^'^''- At each iteration the su- 
perblock consists of a system block, an environment block 
and two extra sites. The initial system and environment 
blocks consist of two sites. The system and environment 
blocks are increased by two sites at a time until a su- 
perblock of = 50 sites is reached. We retain m = 230 
density matrix eigenstates in the basis truncation proce- 
dure. 

The system, environment and super block Hamiltoni- 
ans and density matrices are block diagonalised using 
the total particle number N and total z-spin opera- 
tors. States from different symmetry sectors are found 
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by projecting trial states from the iterative, sparse ma- 
trix diagonalisation procedure into the correct symmetry 
sector by means Q^ithe operators C2h, J and the spin- 
parity operator PEJl. Because the density matrix com- 
mutes with the block J and P operators, the superblock 
states calculated are exact eigenstates of J and P at all 
stages of the calculation. They are also exact eigenstates 
of C2h by construction because the environment block is 
the reflection of the system block. We check the DMRG 
program by ensuring that it reproduces the results of the 
exact diagonalisation program for the first two iterations. 
The exact diagonalisation program itself reproduces the 
N = 8 site results fromQ. 



III. RESULTS AND DISCUSSION 
A. Accuracy of the DMRG calculations 

We begin with the non-interacting limit {U = V ^ 0), 
where we have performed separate exact diagonalisations 
for long chains in order to evaluate the DMRG accuracy. 
The non-interacting case is of particular interest since it 
gives a worst case accuracy. Errors are expected to be 
smaller in the interacting case where particles are more 
localised in position space, the DMRG becoming exact 
in the atomic {t = 0) limit. 

In the DMRG calculations the initial two systems 
{N = 6, 10) are treated exactly. The errors for these 
systems are merely a result of the limited precision of 
the sparse matrix diagonalisation algorithm (the accu- 
racy could be increased by running the programs at 
higher precision). The DMRG truncation error sets in 
at = 14 where the relative error is ^ 10^^. The error 
increases with system size and is ^ 10~^ for N = 30-50. 

In addition to excitation energies, it is important 
to look at quantities such as transition moments and 
particle-hole separations, which pertain to the wavefunc- 
tion rather than the energy. The truncation error for 
transition moments is ~ 10~^ for N = 30-50, i.e. one 
order of magnitude larger than for excitation energies, 
but well within what is acceptable for our purposes. The 
particle-hole separations have similar errors to the tran- 
sition moments. 

The DMRG is a truncated basis method with system- 
atically reducible error. In the interacting case it is im- 
portant to test the accuracy of calculations by varying the 
single source of error, the truncation parameter to, and 
checking for convergence. We asses the error by running 
the program for a number of values of m. Convergence 
results for a number of quantities are given in Table | for 
the TV = 50 site system. It is found that the ground state 
energy converges to 5 or 6 figures and gaps are resolved 
to within 0.01% for the charge gap and IB^ states, with 
slower convergence for the higher excitations where errors 
range up to one or two percent. The errors in transition 
dipole moments are larger but are still small enough to 



make a clear identification of the essential NLO states. 
Errors range from around 0.5% for the lAg — > IB^ tran- 
sition to 1-2% for the mAg IBu transition. The errors 
in the particle-hole separations are found to range from 
0.1% for the lAg state to 1-3% for the higher excitations. 
Convergence is sufficiently good to allow the clear clas- 
sification of states as bound or unbound on the basis of 
the scaling of the particle-hole separation. 

TABLE I. The ground state energy (lAg), the charge gap 
Eg, the IBu and AAg energy gaps, various transition mo- 
ments {jBu\fi\lAg) and particle-hole separations {1)^q for the 
= 50 site system calculated using the DMRG for a number 
of truncation parameter values m. 



m 




Energy 




{jBu\fi\mAg) 


(Oso 


lAg 




1B„ 


4Ag 


IBu 3Bu 


2Bu 


4B„ 


64 


-102.97649 


1.2190 


1.1122 


1.1761 


1.479 9.472 


1.883 


1.942 


100 


-102.97925 


1.2212 


1.1058 


1.5455 


7.151 6.599 


2.018 


1.726 


150 


-102.97969 


1.2213 


1.1044 


1.4801 


7.065 6.847 


2.061 


1.632 


185 


-102.97996 


1.2218 


1.1026 


1.3855 


6.962 7.297 


2.078 


1.602 


230 


-102.98002 


1.2218 


1.1025 


1.3431 


6.945 7.330 


2.085 


1.596 



B. Excitation energies and transition moments 

The evolution of the energy spectra with system size 
is shown in Fig. || and Fig. |j for the Ag and Bu sectors 
respectively. The charge gap is also shown for compar- 
ison. There is exactly one excitation below the charge 
gap for all system sizes in the Bu symmetry sector, and 
for most systems in the Ag sector. The 3Ag drops below 
the charge gap for N = 42-50. 



CP 

u 

0) 1-5 

a 




10 20 30 40 50 



N > 

FIG. 3. Excitation energies (relative to the lAg) in the Ag 
symmetry sector as a function of TV. The charge gap (dashed 
line) is shown for comparison. 

Transition moments between various Ag states and the 
IBu are given in Table ^ The mAg can be identified as 
the 5Ag for N — 10-50. The transition moment with the 
2Ag increases with system size, while transition moments 
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with the 3Ag and the AAg decrease for A'' > 22. The 
charge gap is always below the mAg energy so the charge 
gap criterion will give lower exciton binding energies than 
the mAg criterion. 

TABLE II. Transition moments {jAg\fl\lBu} for various 
states jAg and system sizes A''. 

{jAg\f,\lB^) 



iV 


lAg 


2Ag 


3Ag 


4Ag 


5Ag 


6 


1.529 


0.239 


0.185 


2.697 


0.195 


10 


2.144 


0.664 


0.000 


0.970 


4.006 


14 


2.613 


1.210 


0.179 


2.089 


4.879 


18 


2.987 


1.794 


0.270 


2.592 


5.691 


22 


3.300 


2.369 


0.241 


2.779 


6.384 


26 


3.548 


2.924 


0.173 


2.770 


6.920 


30 


3.789 


3.408 


0.112 


2.639 


7.290 


34 


4.071 


3.706 


0.069 


2.421 


7.464 


38 


4.303 


3.989 


0.046 


2.162 


7.543 


42 


4.528 


4.203 


0.043 


1.526 


7.627 


46 


4.636 


4.703 


0.046 


1.786 


7.594 


50 


4.831 


4.899 


0.058 


1.603 


7.521 
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FIG. 4. Excitation energies (relative to the lAg) in the Bu 
symmetry sector as a function of A'^. The charge gap (dashed 
line) and is shown for comparison. 



argument, based on th e scal ing of the particle- hole sepa- 
ration, is given in Sec. IIIC, 



The transition moments for the identification of the 

There are two qualitative 



nBu are given in Table III 



differences between these transition moments and those 
used for the identification of the mAg (Table ||) . Firstly, 
while the state number of the mAg is basically constant 
(to = 5, apart from N = 6), the state number of the nBu 
changes with the system size when the definition of Sec. 
I is apphed. It is the 5B„ at N=6, the 4B„ at iV = 10, 14 
and the 2B„ at iV = 18, 22, 26. 

TABLE III. Transition moments {jBu\fi\mAg) for various 
states jBu and system sizes A'^. 









{jBu 


fj.\mAg) 






iV 


IBu 


2Bu 


3Bu 


4Bu 


5B„ 


6Bu 


6 


2.697 


0.908 


0.281 


1.916 


2.135 


0.041 


10 


4.006 


2.059 


0.391 


3.977 


0.235 


1.735 


14 


4.879 


3.524 


0.016 


4.640 


0.287 


2.446 


18 


5.691 


5.278 


0.845 


4.771 


3.299 


1.291 


22 


6.384 


6.892 


2.664 


0.396 


5.997 


2.207 


26 


6.920 


6.860 


6.585 


1.822 


0.987 


6.361 



C. Particle- hole separation 



I, 
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FIG. 5. Particle-hole separation (Z)^ of the Ag states as a 
function of A^. The plot shows all states up to the mAg. 



Secondly, for many system sizes, the nBu is ill-defined 
in the sense that there are several transition moments 
close to the maximum value and to that of the nB^- For 
instance, at = 18, the li?„, 2i?„ and 4i?„ have very 
similar transition moments; at A^ = 22, the IBu, 2Bu 
and 5i?„; and at A^ = 26, the li?„, 2Bu, 3i?„, and 6i?„. 
It was found inQ that conduction-band states had strong 
dipole coupling to several neighbouring conduction-band 
states in the dipole coupled symmetry sector. This was 
seen to distinguish band-to-band transitions from transi- 
tions involving excitons. If the mAg is the onset of the 
conduction-band, i.e. the lowest lying unbound Ag state, 
then our results imply that the IB^ is the onset of the 
conduction-band in the B^ symmetry sector. A stronger 



The particle-hole separation (Z)^^ of the Ag states is 
plotted as a function of A^ in Fig. ^. We see that the 
mAg displays completely different behaviour from all the 
lower lying Ag states. For the lower states, the average 
separation converges rapidly to the grouni-state (bulk 
limit) value of around 1.3 chemical bondsESI. By con- 
trast, for the mAg the separation is significantly larger 
for all A^, exhibiting a far stronger A'^ dependence, and 
converging to a different limit from the bulk limit if con- 
verging at all. It follows that the Ag states below the 
mAg are bound and the mAg is unbound. 

The striking difference between the mAg and the lower 
Ag states is also borne out in the behaviour of the reduced 

particle-hole separation {1)%'^^'' depicted in Fig. |[ The 
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FIG. 6. The reduced particle-hole separation (Ojv° of the FIG. 7. Particle-hole separation (l)j^ of the Bu states as a 

Ag states as a function of A'^. The plot shows all states up to function of N. The lAg state is plotted for reference. 



reduced separation of the mAg scales linearly with N, i.e. 
as for completely independent particles and holes. For 
N = 50, the average separation is around 25 chemical 
bonds. The lower states, on the other hand, all have 



(red) 
N 



bounded separation. For instance, for the 2Ag (1) 
scales linearly up to = 18 where it levels off at around 
8 chemical bonds. The interpretation of this is clear: 
the 2Ag is comprised of excitons with an average size of 
8 bonds. When the system reaches a sufficiently large 
size the separation settles at the intrinsic average size of 
the exciton. For smaller systems, around or below the 
exciton size, the average separation is dictated by the 
system confinement. 

In summary, the Ag states below the mAg are bound, 
and the mAg is the lowest lying free charge-carrier state 
in this symmetry sector. We have only performed calcu- 
lations for up to A'^ = 50, and thus in principle we can- 
not be certain that the reduced separation of the mAg 
does not eventually taper off. However, it is clear from 
our calculations that the mAg behaves as a free charge- 
carrier state for systems of up to 50 sites. If an up- 
per bound does exist, then the mAg must be extremely 
weakly bound. The mAg energy is therefore a good ref- 
erence energy when calculating binding energies in this 
symmetry sector. The charge gap, on the other hand, 
falls amongst the bound states. Some of the bound states 
are above the charge gap, while others are below. Thus, 
the charge gap fails to discriminate between bound and 
unbound states. 

We now turn to the dipole-allowed (Bu) symmetry sec- 
tor. The particle-hole separation is shown in Fig. 0. The 
separations depend strongly on N, indicating that there 
are no bound states in this sector. This is confirmed 
by the reduced particle-hole separation, shown in Fig. 
The IBu scales almost perfectly linearly with N and the 
2Bu increases at least as rapidly over the plotted range. 
This explains the absen ce of a well defined nB^ state 
above the (see Sec. IIIB). That is, the IBu is the 
lowest unbound state in the Bu sector. We note that 



the charge gap lies above the incorrectly suggest- 

ing that the IBu is a bound state (see Fig. ^). As was 
the case for the Ag sector, the charge gap is not a useful 
criterion for determining whether states are bound and 
unbound in this symmetry sector. 




N 



FIG. 8. The reduced particle-hole separation (/)^°'*' of the 
Bu states as a function of A''. 



D. lonicity 

The ionicity gives additional insight into the nature 
of the low-lying states. It is instructive to first study a 
case of relatively strong Coulomb interaction {U = 10, 
V = A). As seen in Fig. ^ there is a sudden jump in the 

ionicity at the mAg, for which ^/^ = 1.4. The ground 

state has an ionicity of 0.45, and an excitation to the 
mAg creates almost exactly one particle-hole pair. Lower 
lying states, on the other hand, follow a more continuous 
evolution with similar ionicity to that of the ground state. 
In fact, the ionicity decreases, the {m—l)Ag being almost 
purely covalent. By contrast, there is no such jump in 
the Bu symmetry sector at the n_B„. Instead all states 
have about the same ionicity, similar to that of the mAg . 
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TABLE IV. The average ionic separation of jAg and j_B„ 
states for iV = 10. 



FIG. 9. The ionicity of the jAg and jBu states for U = 10, 
1/ = 4, TV = 10. 



For the weaker, more realistic, interaction (see Sec. 

Although there 



[I A), the ionicity is shown in Fig. |l£ 
are large quantitative differences, many qualitative fea- 
tures remain. In the Ag sector, the ionicity decreases 
continuously up to the mAg where a distinct jump oc- 
curs. Again, the ionicity of the _B„ states remains con- 
stant at around the same value as the mAg. There is, 
however, one qualitative change, namely that the spec- 
trum above the uiAg is now a mixture of states with high 
and low ionicity. A picture where the mAg is the onset 
of a band of states which are all unbound appears to be 
too simplistic, at least from the = 10 site data. 




8 10 12 14 IS 18 20 



FIG. 10. The ionicity of the jAg and jBu states for (7 = 3, 
V = 1.2, N = 10. 



The average ionic separation, shown in Table IV, also 
shows a jump at the mAg, from slightly more than 2 for 
the lower states to 3.5. In the other symmetry sector, the 
maximum ionic separation is 3.8 for the nBu- However, 
the jump is less pronounced since the lower states all have 
relatively large separations of around 3.0. 

To conclude, the mAg excitation involves the creation 
of more particle-hole pairs than the lower lying states in 
the Ag sector. Furthermore, the particle-hole pairs as less 
strongly bound. The large charge separation of the mAg 
is a result of two factors: a redistribution of weights in 
the wavefunction, from covalent to ionic configurations. 











j 


A, 




Bu 


1 


2.23 




2.97 


2 


2.21 




3.21 


3 


2.04 




2.82 


4 = n 


2.11 




3.77 


5 = m 


3.54 




3.08 



and an increased charge separation within the ionic part 
of the wave function. By contrast, in the symme- 
try sector the difference between the nBu and the lower 
states is far less pronoun ced, a result which is consistent 
with the findings of Sec. IIIC. 



IV. CONCLUSIONS 

In this paper we have studied the particle-hole separa- 
tion in the fundamental model of conjugated polymers — 
the dimerised, extended Hubbard model — in order to 
critically assess criteria commonly used to determine 
whether excitations are excitonic or consist of uncorre- 
lated particle-hole pairs. The chosen parameter values 
are typical values used to describe polymers such as poly- 
acetylene within the model. The model was solved for 
the charge gap and a number of states in the ground 
state {Ag) and dipole allowed [By) symmetry sectors us- 
ing the Density Matrix Renormalisation Group (DMRG) 
method for systems of up to iV = 50 sites. It was shown 
that the DMRG can be used to accurately compute en- 
ergy gaps, transition dipole moments and particle-hole 
separations for these states, with relative errors ranging 
from a fraction of a percent for the lowest states to a 
few percent for higher states. The particle-hole separa- 
tions of bound and unbound states were shown to scale 
differently with system size, a fact which can be used to 
discriminate between the states. 

It was found that the charge gap, often used to define 
the binding energy of excitons, is not a useful criterion 
by which to decide whether a state is bound or unbound. 
The results for the scaling of the particle-hole separation 
show that there is an unbound state below the charge 
gap in the Bu symmetry sector, and bound states states 
above the charge gap in the Ag sector. Thus, it is possible 
to have states below the charge gap which are unbound, 
and states above the charge gap which arc bound. In 
fact, the results were in keeping with the picture that the 
essential non-linear optical state mAg marks the onset of 
unbound states in the Ag symmetry sector. On the other 
hand, in the Bu sector it was found that all states were 
unbound for the parameters considered, and that there 
was no well defined nB^ state. That is, the \Bu state 
marks the onset of unbound excitations in the Bu sector. 

Finally, we note that, although the dimerised, ex- 
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tended Hubbard model is a fundamental model which is 
believed to provide a qualitative description of polymers 
such as polyacetylene, one needs to investigate the effects 
of long range Coulomb interactions, interchain coupling 
and electron-phonon effects before one can begin to make 
reliable statements about the nature of excitons in spe- 
cific, real conjugated systems. We are currently develop- 
ing efficient, vectorised DMRG codes with which we will 
be able to study these effects. 
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APPENDIX A: SCALING OF THE 
CORRELATION FUNCTION AND 
PARTICLE-HOLE SEPARATION IN THE 
NON-INTERACTING LIMIT 

In this appendix we present some results for the cor- 
relation functions of excited states of (||) in the non- 
interacting limit {U — V — 0). In arriving at these 
results, we make use of well known analytical results for 
open and periodic systems and diagonalisations of open 
systems of up to 4098 sites. 

The bulk correlation function is given by 



Coo(0 



lim Cn(1) 




1 = 

I > 0, even 
I odd, 



(Al) 



(A2) 



where 



F{l,S,±)^- 



7r/2 



-7r/2 



COS ±iS sin 6 
cos 9 ^ i6 sin 9 



d9 



(A3) 



For < (5 < 1, the model has a gap A = A6t and expo- 
nentially decaying correlations 



as / ^ oo (Z odd), 



(A4) 



1 



2tanh~^5 



(A5) 



For the ground state of a system of size N we have 



C]v'°^(0 = Coo(0 = for aU I even. 



and for odd Z, C^^''(/) approaches Cao{l) exponentially 
fast viz 

ci,°'^(0=C^oo(/)[l + O(e-"^)] {I odd), (A7) 

for some positive constant a which is uniform over I. 
It follows that the ground state particle-hole separation 
converges rapidly, according to where 



a — lim 



e::i(2/-i)|Coo(2/-i) 

e;:iIc^(2;-i)| 



(A8) 



is the particle-hol e sep aration defined from the bulk cor- 
relation function (Al). This can be seen in Fig. |^ where 
(Z)^ is plotted for the lAg state. 




with correlation length 



2050 N = 4098 

FIG. 11. log |Cjv(OI for the IBu state as a function of I for 
iV = 6, 10, 18, 34, 66, 130, 258, 514, 1026, 2050 and 4098 
in the non- interacting case (U = V = 0). Also plotted is 
log |Coo(Oli the bulk (ground state) value (dashed line). Only 
the odd values of I are used. The inset shows the same plot 
on a smaller horizontal scale. 

For the excited states the finite lattice correlation func- 
tion Cn{1) approaches ilie infinite lattice value Coo(0 ^.s 
N ^ oo for any fixed lc3. However, as can be seen from 
Fig. |ll|, where we plot log [Cat (01 for the li3„ states as 
a function of I for various values of TV, the convergence 
is very slow. This is shown in Fig. |l2|, where we plot the 
reduced correlation function C^'^'^\l), defined by (^), as 
a function of for a number of values of I. We see 
clearly that the scaling behaviour of the excited state 
correlation function is 
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Cn{1) - Coo[l) + + ^ 



(A9) 



where c!i^^\l) is the Hmiting, reduced correlation func- 
tion viz 



Cg'^'^'Hl) = lim C 



N^oo 



(rod) 
N 



{I). 



(AlO) 
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^ f(l) 
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ooooooooo 
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FIG. 12. Reduced correlation function C]y° (l) for the l_Bu 
state as a function of 1/A'^ for various (a) odd and (b) even 
values of I in the non-interacting case {U = V = 0). 



Using a polynomial fit of the form (A9) and extrap- 



olating to the N = oo limit, C'^'^\l) is calculated and 



plotted in Fig. 13. We see that Cio (0 decays expo- 
nentially with Z, and vanishes for even I. In fact, the 
scaling form (A9) can be derived explicitly for the case 
of periodic boundary conditions, where it can be shown 
that 



where the functions /(/) and g{l) rapidly approach non- 
zero constants as I is increased (i.e. c!^'^\l) decays expo- 
nentially and has correlation length 2^). The functions 
/ and g, derived from the polynomial fit to (|A^), are 
plotted in the inset of Fig. ^ 

Now, the straightness of the curves in Fig. |lj indi- 
cates that the higher order terms in ( [A9| ) are small, and 
Fig. ^ indicates that (A£) holds for a range of I values 



FIG. 13. Limiting reduced correlation function c!^'^\l) 
for the IBu state as a function of I in the non-interacting 
case {U = V = 0). The inset shows the scaling functions 
f{l) and g{l) {f{l) is only defined for odd values of I because 
C!i'"^\l) = for even 0- 



that is of 0{N). These results, together with the defi- 
nition (H), imply that the particle-hole separation for an 
excited state scales according to (||), with limiting value 
b = limjv^oo (Ojv- This can be seen from Fig. 0, where 
is plotted as a function of 1/iV for the IB^ and 2Ag 
states. We note that b ^ a i.e. the ground and excited 
states have different particle-hole separations in the bulk 
limit, even though their correlation functions approach 
the same limiting fun ction. This is due to the fact that 
the third term in ( |A9| ) docs not decay with I. 
Finally, we note that, from 



( |A9| ) and ( |All| ), the 

scaling behaviour of the reduced correlation function is 



(A12) 



It follows that the reduced separation (1)%'^^'' scales lin- 
early with TV. This is clearly shown in Fig. where 

, (red) 



(All) {1)n plotted for the IBu and 2Ag states. 
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